Influence of the impellers' geometry on turbulent von Karman swirling flows 
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We report numerical evidences of the influence of the blades' shape and of their direction of 
rotation on turbulent flows generated by counter-rotating bladed-disks enclosed in a cylinder, known 
as the von Karman swirling flow. We show that, although the mean flows generated by stirrers with 
straight or curved blades can be topologically similar, turbulence production occurs in different 
zones. The structure and intensity of the turbulent fluctuating field also varies with the geometry, 
and retains significant anisotropy. From the direct numerical simulations fields, we analyse the 
dynamics of enstrophy and we show that the distribution of relative helicity h = u ■ h/(\u\\w\) 
depends on the forcing mecanism of turbulence. 



Introduction - In order to experimentally study statis- 
tically stationary turbulence, scientists developed in the 
last two decades the so-called von Karman device [1-4], 
designed to produce inertial turbulent flows by the rota- 
tion of two facing bladed impellers enclosed in a cylin- 
drical cell. Different regimes for the turbulent flow are 
observed when the rotation rates /i and f% of the two im- 
pellers are varied, and for different aspect ratios Y = H/R 
of the fluid cell, where H is the distance between the two 
disks and R the inner radius of the cylindrical container. 
The kinematic viscosity of the fluid v is also a parame- 
ter, and one builds a Reynolds number Rex = 2irR 2 fv^ 1 
based on the intensity of the forcing by the bladed disks 
/ = (ff + /I) 1 / 2 . However, the variability of the differ- 
ent experimental flows observed can only be explained by 
considering the precise geometry of the blades mounted 
on the rotating disks — for example, a dynamo action 
can be produced with a certain kind of impellers [5] and 
not with other ones [6] for an equivalent Re^f. Many fea- 
tures of the turbulent flows observed with curved blades 
cannot be observed with straight ones: for instance, the 
multistability of the mean flow (between one- or two-cell 
states) is observed with curved blades but nonexistent 
with straight blades [ ]. Moreover, the direction of ro- 
tation of the blades, when curved, actually modify the 
properties of the flows [8, 9]. In addition, the homogene- 
ity and isotropy of turbulence in the core of the cell, once 
thought established, is now questioned. In this letter, we 
analyze the mechanisms of turbulence production from 
the forcing by either straight or curved blades, the zones 
where energy is injected and transported, and how it in- 
fluences the homogeneity and isotropy properties. Such 
an analysis is permitted by direct numerical simulations 
(DNS) with a volume-penalization technique for the ro- 
tating blades and the cylindrical container. We thus have 
access to the phenomena in a way not accessible by ex- 
perimental realizations. 

Flow geometry, parameters and algorithm - To allow 
comparison with existing experiments of von Karman 
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swirling flows, we consider a cylindrical container in 
which two facing rotors are placed, with an aspect ratio 
of the fluid cell T — 1.8, equivalent to [4, 7-9]. The radius 
Rd of the disks is the same as the TM60 impellers used 
in these experimental works (Rd = 0.925i?) as well as 
the height of the blades (Hi, = 0.2R) and their curvature 
radius (R CU rv = 0.5i?) for the curved ones. The main 
difference with the experiments is a number of blades 
reduced to four (6 straight ones in [2], 16 curved ones 
in [7] for example) in the simulations, due to numerical 
constraints, as shown on figure 1. We perform DNS of 
the Navier-Stokes equations using a 3D Fourier pseudo- 
spectral method in a periodic domain. The above ge- 
ometry of the flow and the no-slip boundary conditions 
are imposed by a volume-penalization method [ ] , which 
allows an accurate representation of the fixed cylindrical 
container and of the mobile impellers. Technical details 
and validation of this method are explained in [11]. The 
penalized Navier-Stokes equations solved are therefore 

du 

— - u x w = -Vpt +v\J 2 u - ri~ 1 x(u- u s ) , (1) 
at 

where u and oj — V x u are the velocity and vorticity 
fields, pt is the total pressure and \ is the mask func- 
tion separating the spatial domain into solid (% = 1) and 
fluid (x — 0) regions. The last term on the right-hand- 
side of equation (1) is the volume penalization term and 
the combined specification of x( x ,t) and of the bound- 
ary velocity u s (x,t) sets the required boundary condi- 
tions. r\ is the penalization parameter chosen as small 
as possible (typically 77 ~ 10~ 3 ). In the present work, 
we focus on the exact counter-rotating regime /1 = — /2, 
where / = 0.1 with a fluid viscosity v = 2 x 10~ 4 (all 
the parameters are non-dimensionalized). The radius of 
the cylinder is R = 0.9L/j, where = 2ir is the hori- 
zontal size of the total resolution domain, whose vertical 
extent is L v = An. The Reynolds number is therefore 
Re K « 25000, above the value 20000 at which fully de- 
vcloppcd turbulence is expected [ ], but well below that 
for a bi-stable behavior, at Re^ « 10 5 [12]. The number 
of collocation points in the horizontal and vertical direc- 
tions are N x y = 512 and N z — 1024. We consider three 
different cases: case S with straight blades (left figure 1) 
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FIG. 1. (Color online) Geometries of the simulated blades, 
corresponding to the initial distribution of the mask function 
X in a horizontal plane containing the blades (red=solid; yel- 
low=fluid). For the curved blades, the direction of rotation is 
either negative (concave face forward) or positive © (convex 
face forward). 
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FIG. 2. (Color online) Mean flows in a half vertical plane 
(symmetry axis on the left, solid wall on the right): C- for 
curved blades in negative direction, C+ for the positive di- 
rection, S for straight blades. In color the azimuthal velocity 
component ug ; the vectors are projections of the mean veloc- 
ity field on the rz-plane. Re^ ~ 25000 for all cases. and 
show the rotation direction of the top and bottom impellers, 
located in the white spaces. 



and case C- (resp. C+) with curved blades rotating in 
the negative (resp. positive) direction of rotation (right 
figure 1). The stirrers' rotation is progressively increased 
from rest to avoid a discontinuity in time. 

Statistically stationary states - Although at the same 
Reif, the different geometries induce important differ- 
ences in the flow evolution and in the eventual statisti- 
cally stationary state. First, the mean total kinetic en- 
ergy in the fluid cell is k — (it 2 ) /2 ~ 0.1 for case S, 
but reduced to k ~ 0.06 for case C+ and to k ~ 0.03 
for case C-. Consequently, the energy cascade from the 
large forced scales to the dissipative ones is reduced, so 
that the steady state is reached very shortly in case 
C- {t s tatf — 0.5), but over a longer time in cases S 
(t stat f ~ 3.5) and C+ (t sta tf ^ 6). Figure 2 shows 
the corresponding computed mean flows, obtained by az- 
imuthal and temporal averaging from t — t stat to the 
end of the computations. Clearly, the three mean flows 




FIG. 3. (Color online) Volume rendering of the instantaneous 
enstrophy field Z — uj 2 /2 for the straight blades geometry 
(left) and for the positive-rotating curved blades one (right). 
The full spatial numerical domain is also delineated. 



are significantly different, especially for case C- where 
the curved blades rotating in this reverse direction are 
not able to generate the same large pair of toroidal vor- 
tices in the fluid domain, as cases S and C+, which agree 
with experimental observations for a symmetric forcing 
fx = — / 2 . Accordingly, case C- has too low energy to be 
considered turbulent — the rms velocity in the core re- 
gion is (( u2 ) core ) 1 ^ 2 ~ 10 -4 , compared to the maximum 
velocity of the impellers it s , m ax = 27r/i? w 1.8 — and 
is not discussed any more from here on. The two cases 
S and C+ exhibit noticeable differences: the magnitude 
of azimuthal velocity component is much larger for the 
straight blades case S than for the curved one C+. This 
dominance is however counter-balanced by larger verti- 
cal velocities towards the disks in the neighborhood of 
the impellers at radii r* = r/R < 0.5 for case C+. In 
the curved blades simulation, one also notices a thin re- 
gion r* > 0.85 which does not belong to the recirculation 
cells. In this region, when r* increases, the azimuthal 
velocity component is depleted slightly before reaching a 
local maximum and connecting with the inner boundary 
layer of the cylinder (wall at r* = 1). The C+ and S 
toroidal cells therefore have a different shape and size, 
and a different energy distribution. 

Turbulence properties - In the core region, the flow 
corresponds to inertial turbulence. The Reynolds num- 
ber based on the Taylor micro-scale is Re>, ~ 120 for the 
straight blades (case S) and ReA ~ 80 for the curved 
blades (case C+), although Re^ is unchanged. This 
demonstrates that the turbulent field intensity depends 
importantly on the geometry of the blades at equivalent 
rotation frequency. The small scale turbulent structures 
are depicted on figure 3 with isosurfaces of the enstrophy 
Z = uj 2 /2 for cases S and C+. These figures from DNS 
fields allow an original characterization of unexpected 
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turbulence features, not accessible to the experimental 
observations: the intense turbulent vortices are clustered 
around the horizontal symmetry plane at z = for the 
straight blades simulation, showing that most of the vor- 
ticity production is due to the horizontal shear layer of 
azimuthal velocity; for the curved blades flow, the distri- 
bution of enstrophy highlights intense vortical structures 
along the axis of rotation and the horizontal clustering 
is almost not observed in this configuration. The mech- 
anism of turbulence generation seems to vary strongly 
with the shape of the blades. We confirm this by consid- 
ering the dynamical equation of enstrophy Z 

T)tZ = uiiLOjSij + vu}i\7 2 0Ji, (2) 

where Sij — (diUj + djUi)/2 is the rate of strain tensor 
and D t = d t +u- V. The first term on the right-hand-side 
of equation (2) is the vortex-stretching term and is prin- 
cipally responsible for enstrophy production. The second 
term linked to viscosity can also be a source of vorticity 
and enstrophy production by the tilting of vortices or by 
helping reconnection (see [13]), but this is globally not 
the case for these flows, as can be seen on figure 4 (very 
few occurences of i/UiW 2 uj i > 0). This figure shows the 
vortex stretching term and the viscous term for cases S 
and C+: we first observe that the straight blades pro- 
duce highly turbulent small structures close to their ex- 
tremities (r* — > 1), even when the stationary state is 
reached. This mechanism is different in run C+ since it 
is no longer observed once the mean recirculation cells 
are established. Second, with straight blades, enstrophy 
production is localized in the horizontal shear layer and 
the fluctuations are advected by the mean flow towards 
the center of the fluid domain. This phenomenon changes 
with curved blades, where the vortices are created by the 
large negative radial velocities close to the center and by 
the vertical shear induced by azimuthal velocities of op- 
posite sign within the recirculation cells. This explains 
why the turbulent structures are almost confined to a re- 
gion along the vertical axis of rotation. The generation of 
turbulence is thus very different according to the shape of 
the blades, which also explains the different total energy 
and mean flow topology. 

Anisotropy of turbulence - We now refine the char- 
acterization of turbulence in the center of the fluid do- 
main. We define a core region T> c where r* < 0.3 and 
2\z\/H < 0.3 and we compute from the S and C+ veloc- 
ity fields the Taylor micro-scale in each direction defined 
by 

A,- ((4fc)/(3(d 4Ul ) 2 )) 1/2 (3) 

(no summation over i). We previously observed that the 
Reynolds number Re^ is larger in the straight blades case, 
but this is essentially due to a larger kinetic energy, the 
Taylor scales A.; being lower (A z ps 0.38 for the case S 
and A z ps 0.42 for the case C+ for example). In both 
cases, the vertical length scale A z is greater than the hor- 
izontal ones X x .y The difference is about 5% in case S 




FIG. 4. (Color online) Instantaneous distribution of the 
vortex stretching term LOitOjSij (top) and the viscous term 
voOiX7 2 u!i (bottom) of the evolution equation of the enstrophy 
(2) for the straight blades (left) and the curved blades (right). 
Both quantities are non-dimensionalized by (Z) (s) /k, where 
(e) is the mean dissipation rate. 



and reaches about 15% in case C+. The separation of 
the directional Taylor scales shows that the flow is not 
isotropic in the small scales, also confirmed by the plot 
of the two-point vertical vorticity correlation function 

(uj z (x)uj z (x + lei))^ 
K?(l) = , Pc , (4) 

\j{ul{x)) Vc {u) 2 z {x + lei)) Vc 

where < • >t> c denotes spatial averaging over D c . The 
vorticity correlation function differs from the velocity 
one, since in the latter the mean flow contribution, al- 
though relatively weak, prevents to estimate accurate 
turbulent correlation lengths (one could also subtract this 
mean flow in the velocity field, although we have not done 
so here). Figure 5 shows that vorticity is much more 
correlated in the vertical direction than in the horizon- 
tal ones, which is consistent with the computed Taylor 
scales, the vertical ones being larger than the horizon- 
tal ones. It also clearly emphasizes the anisotropy of 
turbulence in the central region, for all blades geometry, 
contrarily to some initial assumptions when experimental 
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FIG. 5. (Color online) Correlation function based on the ver- 
tical vorticity component to z for a horizontal or a vertical 
separation, of the flows generated by straight (S) or curved 
blades (C+). 
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FIG. 6. (Color online) Pdfs of cos# — where 9 is the angle 
between velocity and vorticity vectors — equal to the 'relative 
helicity' h. The solid black line denotes the results of isotropic 
simulations by Rogers & Moin [16] and the dotted grey one 
those of Kerr [15] with a large-scale forcing. 

characterization of von Karman swirling flows began. 

Beltramization - We consider the relative orientation 
of velocity and vorticity (see [14, 15]) by computing the 



relative helicity h = u-u>/\u\\u>\ = cos(9, 9 being the 
angle between the two vectors. This angle is also im- 
portant because it measures somehow the amplitude of 
the non- linear term in equation (1) and the efficiency of 
the energy cascade. We compute the probability density 
function (pdf) of h and compare it to previous observa- 
tions in simulations of homogeneous isotropic turbulence 
[15, 16]. The results presented on figure 6 show surpris- 
ing features. On the one hand, in the straight blades 
case S, the pdf of h is approximately flat, as if veloc- 
ity and vorticity were uncorrelated, thus differing from 
the homogeneous isotropic turbulence pdf which curves 
up at cos 9 — ±1. On the other hand, run C+ shows 
that the flow enhances velocity-vorticity alignment and 
inhibits anti- alignment. This is explained by the pres- 
ence of the mean flow which is non-negligible even in the 
region T> c , as shown by figure 2 for case C+. As the 
vorticity is mainly vertical (see figure 3 with the intense 
vertical vortices) and more likely to be aligned with the 
rotation of the closer disk, and as the vertical velocity 
is linked to the mean recirculation cells, we believe that 
this alignment is not a feature of turbulence but a trace 
of the forcing by the mean flow. Opposite rotation of the 
disks would lead to the anti-symmetric distribution with 
a peak at cos 9 = —1. But we also note that a similar 
distribution of the pdf of h is also observed by [15] in sim- 
ulations of isotropic turbulence with large-scale forcing. 
Overall, this confirms that forcing turbulence modifies in 
some way its 'natural' dynamics and structure, and that 
high Reynolds number homogeneous isotropic turbulence 
may be very hard to achieve. 

Conclusion - We have highlighted the importance of 
the shape of the blades on the generation of turbulence, 
and on the turbulence itself. Differences are observed 
for the mean flows even at identical Rex, since turbu- 
lence production can appear in different zones of the cell. 
Visualizations, directional Taylor micro-scales and corre- 
lation functions show the non negligeable anisotropy of 
turbulence in the core region. 
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